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Internal  Wave  Generation  by  Submerged  Bodies:  Mean  Flow  Effects 


Introduction 

The  passage  of  a  body  thru  a  density  stratified  fluid  presents  a 
complex  localized  disturbance  to  the  system.  It  is  a  gener&j.  property 
of  such  systems  that  disturbance  energy  is  propagated  away  from  the 
source  region  in  the  form  of  an  internal  gravity  wave  field.  However, 
in  many  circumstances,  the  complexity  of  the  disturbance  created  by  a 
real  body  moving  within  a  real  fluid  effectively  destroys  the  possibil¬ 
ity  of  an  exact  mathematical  description  of  the  radiated  wave  field, 
and  we  must  be  content  with  an  attempt  to  distill  the  essential  physics 
of  the  phenomenon.  That  is,  we  must  isolate  and  understand  those  mech¬ 
anisms  which  characterize  the  disturbance  and  its  associated  wave  field, 
particularly  those  mechanisms  which  transcend  individual  pecularities 
of  particular  bodies  and  environments,  in  the  hope  that  the  physical 
picture  which  emerges  will  provide  a  useful  and  viable  description  of 
the  real  phenomenon. 

On  a  primitive  level,  the  disturbance  may  be  resolved  into  three  < 

categories  on  the  basis  of  temporal  dependence.  In  nature  we  must  deal 
with  real  fluids,  so  that  for  physically  interesting  situations  thf  dis¬ 
turbance  will  be  characterized  by  turbulence  arising  from  boundary  layer 
instability.  The  essence  of  turbulent  motion  is  randomness,  and  all 
aspects  of  the  phenomenon  which  require  a  stochastic  description  will  be 
designated  as  "random  effects".  The  remainder  of  the  motions  associated 
with  the  disturbance  and  the  radiated  wave  field  may  be  designated  as 
steady  or  "mean  field  effects"  or  unsteady  or  "transient  effects".  If 
the  body  is  in  steady  horizontal  motion  over  a  flat  bottom,  then  in  a 
reference  frame  moving  with  the  body,  all  non- random  effects  will  be 
mean  field  effects.  In  such  a  frame  of  reference,  transient  effects 
arise  only  as  a  result  of  unsteadiness  in  the  motion  of  the  body  or 
apparent  unsteadiness  associated  with  the  boundary  configuration. 

The  mean  field  effects  associated  with  a  submerged  body  admit  a 
three-fold  classification.  The  first  effect  is  that  which  would  be  ' 

created  by  the  body  in  the  complete  absence  of  viscous  effects.  The 
disturbance  is  associated  with  the  fact  that  the  body  must  systematically 
displace  and  replace  fluid  as  it  moves,  and  will  be  referred  to  as  the 
"hull  effect"  in  the  sequel.  The  second  effect,  the  so-called  "wake 
collapse  effect",  arises  because,  in  most  physically  interesting  situa¬ 
tions,  a  turbulent  wake  forms  in  the  lee  of  the  body.  Mixing  within  the 
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wake  results  in  an  effective  perturbation  of  the  ambient  density  field, 
which  brings  into  play  relatively  strong  buoyancy  forces  tending  to 
destroy  the  density  perturbation  by  flattening  or  collapsing  the  wake 
region.  The  presence  of  a  turbulent  wake  also  presents  the  possibility 
of  an  effective  "displacement  body"  associated  with  streamline  displace¬ 
ments  outside  of  the  turbulent  region  arising  from  momentum  redistribu¬ 
tion  within  the  wake  region.  This  third  effect,  the  "displacement 
thickness  effect",  is  distinct  from  the  wake  collapse  phenomenon,  being 
similar  in  its  internal  wave  making  characteristics  to  the  hull  effect, 
although  its  interaction  with  the  ambient  density  field  may  result  in 
the  distortion  of  the  turbulent  wake  structure. 

All  of  the  mechanisms  mentioned  above  are  capable  of  affecting  the 
radiated  internal  wave  field  in  some  way,  the  exact  details  of  the  wave 
field  depending  to  some  extent  on  the  exact  details  of  the  body  geometry, 
unsteady  motions  and  the  ambient  stratification.  Because  of  the  uncer¬ 
tainty  involved,  exact  comparisons  of  calculated  waveforms  are  of  limited 
utility;  in  order  to  assess  the  relative  importance  of  the  various  con¬ 
tributing  causes,  their  effects  must  be  characterized  by  some  global 
property  which  is  somewhat  insensitive  to  the  exact  details  of  the  flow 
structure.  For  this  reason,  we  choose  the  rate  of  internal  wave  energy 
production  as  a  convenient  index  of  the  source  strength,  and  explicit 
internal  wave  fields  are  not  computed.  It  is  hoped  that  the  comparison, 
based  on  this  global  property,  of  the  relative  importance  of  the  various 
contributing  mechanisms  as  a  function  of  the  sensible  parameters  which 
characterize  the  system,  will  contribute  to  the  basic  understanding  of 
the  physics  of  the  phenomena  associated  with  the  passage  of  a  body  thru 
a  density  stratified  fluid. 

The  hull  effect 


The  passage  of  a  solid  body  thru  a  fluid  environment  is  accompanied 
by  a  continuous  readjustment  of  fluid  elements.  The  bow  end  continu¬ 
ously  displaces  fluid  while  fluid  must  continuously  refill  the  void  left 
by  the  motion  of  the  tail  end,  and  the  problem  is  then  to  understand  the 
processes  involved  in  the  interaction  of  this  disturbance  with  the  ambi¬ 
ent  density  stratification  of  the  environment.  Most  of  the  basic  physics 
of  this  phenomenon  is  contained  in  the  standard  mountain  lee  wave  problem 
for  which  an  extensive  literature  exists  (see  reviews  by  Queney  et  al, 
1960,  Miles,  1969,  Zeytounian,  1969a,  b,  Vergeiner,  1971,  Lilly,  1972: 
see  also  Bretherton,  1969).  Reducing  the  phenomenon  to  its  pm-'st  form, 
we  consider  the  steady  horizontal  potion  of  a  body  in  an  unbounded  uni¬ 
formly  stratified  Boussinesq  fluid  .  The  problem  may  be  solved  by 

$ 

In  the  Boussinesq  approximation  (see,  for  example,  Phillips,  1966,  sec. 
2.4)  density  variations  are  considered  important  only  insofar  as  they 
determine  the  buoyancy,  which  is  consistent  with  our  program  of  retain¬ 
ing  only  the  essence  of  the  problem. 
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following  Dorodnitsyn  (see  Tareyev,  1965)  and  replacing  the  body  by  an 
equivalent  forcing  function  which  stratifies  the  boundary  conditions. 
Conceptually,  the  simplest  relevant  forcing  is  that  provided  by  a  moving 
dipole,  for  which  solutions  have  been  presented  by  Crapper  (1959), 

Crosch  (1964)  and  Miles  (1971a).  The  general  formalism  for  the  solution 
of  problems  of  general  localized  forcing  in  a  stratified  flow  has  been 
developed  by  Lighthill  (1967)  and  Liu  and  Yeh  (1971).  However,  we  will 
approach  the  problem  in  the  classical  mountain  lee  wav.e  tradition,  devel¬ 
oping  an  integral  transform  solution  to  the  boundary  value  problem*.  As 
noted  previously,  there  is  a  vast  literature  on  the  classical  lee  wave 
problem,  and  therefore  we  will  give  a  rather  concise  derivation  of  the 
solution.  A  thorough  discussion  of  integral  transform  techniques  is 
given  by  Sneddon  (1972). 

Although  the  full  nonlinear  problem  for  a  two  dimensional  object 
such  as  a  ridge  or  horizontal  cylinder  aligned  perpendicular  to  a  uni¬ 
form  stratified  flow  may  be  solved  by  invoking  Long's  model  (Long.  1953, 
1955,  see  also  Miles,  1969,  McIntyre,  1972),  the  full  three  dimensional 
problem  remains  intractable  unless  linearized  by  the  assumption  that 
the  disturbance  is  a  small  perturbation  to  the  basic  state.  A  measure 
of  the  relative  importance  of  nonlinear  effects  is  provided  by  the  in¬ 
verse  internal  Froude  number 

F  ~/  -  (i) 

u 

where  U  is  the  speed  of  the  oncoming  stream  (the  speed  of  the  body  in  a 
stationary  frame  of  reference),  H  is  a  characteristic  value  of  the 
vertical  dimension  of  the  body  (the  maximum  radius  of  a  body  cf  revolu¬ 
tion)  ,  and  N  is  the  Brunt-Vaisala  frequency 


_  d. 

/° 


where  g  is  the  acceleration  of  gravity  and 


L  djO 

/»  J* 


is  the  vertical 


logarithmic  density  gradient  in  the  basic  stat.-*.  Conceptually,  the 
parameter  F*1  is  the  square  root  of  the  ratio  -jf  the  work  required  to 
raise  an  element  of  fluid  to  the  top  of  the  obstacle  to  the  kinetic 
energy  density  of  the  oncoming  fluid.  The  availability  of  the  finite 
amplitude  solution  for  the  two  dimensional  problem  permits  the  evalua¬ 
tion  of  the  validity  of  the  linearized  solution  as  a  function  of  F"1 
for  that  problem  (see  Miles,  1969,  1971b).  These  two-dimensional 
results  indicate  that  the  linearized  problem  should  provide  a  reasonable 


A  somewhat  more  general  approach  to  the  boundary  value  problem  for  an 
isolated  body  in  a  stratified  fluid  is  considered  by  Mackinnon  et  al 
(1969). 
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for  smaller  values  of  F"1  for  certain  body  shapes, 
dimensional  solution  in  the  limit  of  large  F_1  has  been  discussed 


approximation  for  F  ”  less  than  1/4  or  1/2  although,  as  will  be  discussed 
later,  some  error  may  be  introduced  by  the  linearized  boundary  conditions 

The  approximate  three 
by 

Drazin  (1961)  and  Grimshaw  (1969),  and  that  for  very  small  F"-1  by 
Hawthorne  and  Martin  (1955)  and  Drazin  (1961). 


The  linearized  momentum  balance  equations  for  an  inviscid,  Boussin- 
esq  fluid  in  a  reference  frame  moving  with  the  body  (Fig.  1)  are 


(J  d-±<-  r  -  - 

$x  ft,  a  x 


(3a) 


ax 


jL  OP 
ft  ay 


(3b) 


ax 


& 

o>2 


(3c) 


where  U  is  the  speed  of  the  body,  ^  is  a  reference  density,  and  u, 
w,  ,  p  are  perturbations  of  the  velocity  components,  density,  and 
pressure.  Mass  conservation  and  incompressibility  are  expressed 
respectively  by 


v, 


u¥*+»£-° 

d  u  £_y  y  d  w  ,  q  (5) 

d  x  ay  a  ■?  ' 

where^O  is  the  basic  density  profile,  which  is  a  linear  function  of  the 
vertical  coordinate.  Equations  (3)  -  (5)  may  be  combined  to  form  a 
governing  equation  for  the  amplitude  function,  , 


-  o 


(6) 


where  l/oq/ax  =  w  ,  is  the  Laplacian  operator  and 

is  the  horizontal  Laplacian  (  aVsxa+  e>  Vay*) .  Symmetry 
arguments  indicate  that,  if  the  body  is  symmetric  about  a  horizontal 
plane,  the  problem  may  be  restricted  to  a  semi-finite  domain  formed  by 
considering  the  symmetry  plane  as  a  rigid  surface.  Boundary  conditions 


4 


on  n  are  that  Xf  be  Fourier  transformable  in  the  horizontal  directions, 
and  that  ^  (x,y,0)  =  h(x,y),  where  h(x,y)  is  the  half-body  shape  (includ¬ 
ing  the  rigid  surface  on  which  the  half-body  rests).  The  upper  z-boundary 
condition  is  provided  by  a  radiation  condition  of  energy  propagation  away 
from  the  disturbance,  as  required  by  our  notion  of  causality  (see,  for 
example,  Atkinson,  1949). 

To  solve  equation  (6),  we  introduce  the  two  dimensional  Fourier 
transform 


ce 

*/ 


-  i  (XX  +  Xy) 

(  )  &  o/x  d\ 


so  that  equation  (6)  is  transformed  to 


rj"  *  (x2-*  IV  (£bj*-0*i  =  ° 


where  primes  indicate  derivatives  with  respect  to  z.  The  solution  which 
satisfies  the  boundary  condition  at  z  =  0  is  then 

$(*,*,2)  =  e  </',a  (9) 

where  ytf  A  ~  (X*+  Aa)  ~  ()  .  The  vertical  wave  number  is  rendered 

determinate  by  the  radiation  condition.  For  Xat/a>A/?we  must  have 
f)  ->  O  as  2  -*  oo  ,  while  for  X*£/a<  A/*  the  vertical  compon¬ 
ent  of  group  velocity  3  O' / d/4  »  where  q-  X*C/*~  (X*  + 

must  be  positive,  so  that  X*  *  2  a 

/x  --  s3n  (££*  -0  ;  XaO*<A/3 

z.  / >/v"*(io) 

Inverting  the  transform,  we  than  have 

(f  4  t  (MX  -tjyxyxe) 

t(<x,y, ?•>  ■=  U>  (*,»  <?  t/jr c/Jai 

—  OO 

To  the  linear  approximation,  the  horizontal  force  on  the  body 
associated  with  the  wave  field  (the  wave  drag)  is  given  by 


Z>  *  ff p(*>y>0)  V hfXsy)  dx  Jy 
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and  thus  the  rate  of  working  against  the  drag  or  the  rate  of  internal 
wave  energy  production  is  given  by 


oo 


€  =  U'D  -  Ujjp(x,y>°)  dxJy  as) 

By  employing  the  equation=Tof  motion  (3)  -  (5)  and  Parseval’s  relation 
(see,  for  example,  Bremmerman,  1965,  sec.  14.5)  the  expression  (13)  may 
be  transformed  to 


ex* 


./if' 


x 


($-*')*  JxJ* 


(14) 


e  =  J  )*>  (x**AV* 

—CX3  - 

The  expression  (14)  applies  only  to  the  upper  half  of  the  body.  Restrict¬ 
ing  attention  to  symmetric  bodies,  we  then  have  that 

**  -*  'A 


O  p 

for  the  whole  body. 


(l7*~X*)  dx d A 


(15) 


If  we  consider  an  ellipsoidal  body  such  that 


■X  -~y 

then  we  have  that 


h  (x,  ))  -  j  j  j  /_*  -y*  -  X  *Ly  /l  x  J 


(16) 


'A 


coi  xx  cos^y  dyc/x  (17) 

which  may  be  evaluated  (Erdelyi  et  al,  1954,  pp.  11  &  57)  as 

Am ^ <x*+ *VX*>*1 


(18) 


where  we  have  set  R  =  L  /L  ,  the  aspect  ratio,  and  L  =  L  .  J-j  is  the 

x  y  x  v 

Bessel  function  of  the  first  kind  of  order  V  .  The  transform  (18) 

may  be  simplified  somewhat  by*  noting  that  the  Bessel  functions  of  half 
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integral  order  may  be  expressed  in  terms  of  trigonometric  functions, 
that 


[LtM**A*/t I*)''*] 

(H^r/Ra)‘A 


-  cos  Jj-  {X^s)V/?V 


SO 


(19) 


which  may  be  inserted  into  the  energy  production  integral  (15),  yielding 

e-*A(J  **  ItM 


o  o 


-  x*j  j 

}*)**  (x*+A*/A*)*  J 


fs;o[L  (X**iVA*)'AJ  n  r,*L  (20) 

<  /  (X*  V /?*/'*  J 

The  expression  (20)  is  nondimensionalized  by  defining 

/  =  xu/as 

l  = 

"K/Ti  =  ASL/O 

and  introducing  the  normalization  factor  (Miles,  1969a) 

£  =  e/g/g  0*A//AA. 

where  A0  is  the  frontal  area,  equal  to  ^  //Z  //? 
ellipsoid.  In  normalized  form,  then,  we  have 

,,  rr  ,-///  A*  (/-A*)'a 

e  *  ~*  ( 4*)  j i0c7A*+i*)v*(k**ni' 


for  the 


fr¬ 


et  o 


fs/ojpjA^nV 


i* )  '/* 


-cos[%  ik^rr'jTcJkJi 


(21) 


%  2 V 

for  the  rate  of  internal  wave  energy  production  by  the  steady  horizontal 
motion  (time  scale  =  L/U)  of  an  ellipsoid  with  length-width  aspect 
ratio  R  in  a  uniformly  stratified  fluid  (time  scale  =  l/N) . 


The  energy  production  integral  (21)  may  be  evaluated  numerically 
for  various  values  of  R  and  /  7%  .  Fig.  2  is  a  comparison  between 

the  energy  production  rate  for  an  ellipsoid  with  R  =  1  with  that  for 
other  three  dimensional  bodies  for  which  computations  have  appeared  in 
the  literature  (Blumen,  1965,  Macklnnon  et  al,  1969).  A  striking  quali¬ 
tative  difference  between  the  curves  is  apparent  for  '%/'%<  ^At  the 
production  rate  for  the  ellipsoid  being  markedly  larger  than  that  for 
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the  other  bodies  within  this  range,  with  the  curvature  of  the  E  vs 
curves  of  opposite  sign.  The  structure  for  ellipsoids  of  various  P.  is 
illustrated  for  this  range  of  in  Fig.  3.  This  qualitative 

departure  is  tentatively  attributed  to  -he  breakdown  of  the  linearized 
boundary  condition  around  the  lateral  edges  of  t.,e  ellipsoid.  The 
bluff  sides  of  the  ellipsoid  require  infinite  vertical  velocities  at  the 
edge  in  linear  theory,  while  the  other  bodies,  although  of  infinite 
extent,  are  smooth  in  the  appropriate  lateral  directions,  resulting  in 
finite  vertical  velocities  which  tend  uniformly  to  zero  as  the  height 
parameter  tends  to  zero  .  A  qualitative  difference  is  also  apparent  in 
Fourier  space,  where,  because  of  the  slope  discontinuity,  the  lateral 
transform  of  the  ellipsoid  falls  off  as  ^  ,  while  that  for  the 

other  (infinitely  smooth)  bodies  falls  off  exponentially.  This  slow 
decay  associated  with  the  slope  discontinuity  may  in  fact  dominate  the 
behavior  for  small  However,  as  will  become  appar'  later, 

the  behavior  of  the  energy  produc+J  .  \te  curve  in  this  ameter 
range  is  not  of  crucial  importan  '  _le  that  for  larger  V’_  /rt  , 
which  is  illustrated  in  Fig.  4  for  various  values  of  R,  appears  to  be 
reliable,  provided  the  height  scale  H  of  the  body  is  sufficiently 
small  so  that  F"^^  1/2.  For  prolate  spheroids  (ellipsoids  of  revolu¬ 
tion),  this  requirement  reduces  to  7u  ;  R.  Although  there  is 

no  a  priori  justification  for  the  extrapolation  of  two  dimensional 
results  to  the  three  dimensional  case,  it  should  be  noted  that,  in  the 
two  dimensional  problem,  the  effect  of  at  least  mild  nonlinearity,  such 
that  upstream  influence  and  separation  (see,  for  example,  Miles  1971) 
are  not  brought  into  play,  is  to  increase  internal  wave  energy  produc¬ 
tion  rates  above  those  predicted  by  linear  theory  (Miles,  1969).  For 
very  strong  nonlinearity  (  7^  /7J  >>  R)  we  expect  a  reduction  in  the 
energy  supply  to  the  wave  field,  since  in  this  limit  the  flow  is  con¬ 
fined  predominantly  to  horizontal  planes  (Drazin,  1961,  Grimshaw,  1969). 


Wake  collapse  effect 


In  a  real  fluid,  as  the  flow  proceeds  around  the  body,  a  boundary 
layer  of  retarded  flow  develops  which,  in  the  cases  of  physical  interest, 
rapidly  becomes  unstable  and  degenerates  into  turbulence.  Within  this 
turbulent  boundary  layer  and  its  extension  as  a  turbulent  wake,  turbu¬ 
lent  diffusion  acts  to  reduce  the  ambient  density  gradient,  creating  a 
gravitationally  unstable  situation,  which  is  rectified  by  the  release 
of  the  stored  potential  energy  into  kinetic  energy  of  collapse  as  the 
wake  is  flattened  vertically.  Disturbance  energy  associated  with  this 
collapse  process  radiates  into  the  surrounding  fluid  in  the  form  of  an 
internal  gravity  wave  field.  Although  some  theoretical  progress  based 
primarily  on  traditional  dimensional  and  similarity  arguments,  has  been 
made  on  the  problem  of  turbulent  wakes  in  density  stratified  media 


(Onufriyev, 


1971;  see  also  Monin  and  Yaglom,  1971,  ch.  4), 


j|c 

For  the  bodies  considered  by  Mackinnon  et  al  (1969),  the  height 
parameter  is  in  fact  the  half-width  of  the  body.  See  appendix  for  a 


discussion  of  an  approximate  approach  which  removes  the  difficulty 


associated  with  the  bluff  sides  of  the  ellipsoid. 
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a  complete  description  of  the  process  of  wake  generation,  growth,  and 
collapse,  including  the  interaction  of  the  wake  with  the  surrounding 
fluid,  is  not  available.  We  can,  however,  obtain  bounds  on  the  rate  at 
which  energy  is  supplied  to  the  wave  field  by  invoking  rather  basic 
similarity  and  dimensional  arguments  and  relying  rather  heavily  on 
available  empirical  data. 

Wake  collapse  is  a  continuous  process,  active  whenever  buoyancy 
induced  potential  energy  exists  within  the  wake  region-.  However, 
wake  potential  energy  is  created  at  thr  expense  of  turbulence  kinetic 
energy,  which  is  in  turn  fed  by  mean  flow  velocity  gradients  within  the 
wake  region,  so  that  sensible  creation  of  wake  potential  energy  is  con¬ 
fined  to  a  finite  distance  behind  the  body.  This  aspect  of  the  phenom¬ 
enon  is  reflected  in  the  observation  that  a  dyed  wake  expands  to  a 
maximum  vertical  size  and  then  flattens  out.  After  the  passage  of  the 
point  of  maximum  vertical  expansion,  we  can  expect  that  sensible  poten¬ 
tial  energy  production  is  sharply  curtailed,  so  that  the  amount  of 
potential  energy  created  during  the  expansion  stage  should  be  proportional 
to  the  rate  of  internal  wave  energy  production. 


We  may  estimate  the  amount  of  available  potential  energy  by  con¬ 
sidering  the  wake  cross  section  at  the  point  of  maximum  vertical  expan¬ 
sion.  As  a  model,  we  may  consider  an  ellipsoidal  figure,  with  horizontal 
major  axis,  within  which  there  is  a  uniform  stratification  with  Brunt- 
Vaisala  frequency  Ni,  th6  surrounding  fluid  being  characterized  by  the 
Brunt-Vaisala  frequency  N,  so  that  we  may  define  a  mixing  factor 

■f  -  N  Z  ~  Ni  (22) 

In  order  to  account  for  potential  energy  which  ha3  been  lost  during  the 
expansion  stage,  we  replace  the  ellipse  by  a  circle  with  equal  area. 

The  potential  energy  created  per  unit  wake  length  may  then  be  estimated 
as  the  work  required  to  create  this  partially  mixed  circular  area  of 
fluid  by  distorting  the  ambient  density  field.  Referring  to  Fig.  5, 
the  displacement  J  of  a  fluid  parcel  from  its  equilibrium  position 
is  given  by 


The  buoyancy  force  acting  on  a  fluid  element  whose  equilibrium  position 

F  s  -2}  (24) 


is  z  is  given  by 
£ 


so  that  the  potential  energy  density  within  the  wake  cross  section  is 
given  by  2  g  +& 

A.  =  /°  A/-*  (  U  -  Be)  Ja 

(25) 

/S  A'V* 


_/ 
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and  the  total  potential  energy  per  unit  wake  lergth  is 


p  &  /y*  [ 

A  J# 


where  the  integral  is  over  the  cross  sectional  area.  For  the  circular 
partially  mixed  region,  we  thus  have 

r 


P-  j 


where  r  is  the  radius  of  the  circle.  The  energy  per  unit  wake  length 
supplied  to  the  internal  wave  field  should  then  be  proportional  to 
this.  Choosing  the  constant  of  proportionality  to  be  unity,  the  rate  of 
internal  wave  energy  production  by  the  wake  collapse  mechanism  is  approx¬ 
imated  by 

f  s  ZT/g  A/*t/ 

f  (28) 

Normalizing  as  before,  we  then  have  that 

P  =  -L  A/  /  *  rj 

w  */7r  //  A?  /-/3  (29) 

where  H  is  the  body  radius.  If  H  is  a  characteristic  length  scale 
for  the  wake,  such  as  the  body  raSius  (H0  =  H  in  this  case)  or  a  pro¬ 
peller  radius  in  the  case  of  self-propulsion,  then  we  have  that 


qrr 


where  R  =  L/H  and  fE  =  r2/^2 ,  the  ratio  of  the  wake  cross  sectional 
area  at  maximum  vertical  expansion  to  the  characteristic  wake  area, 
order  to  proceed  further,  we  must  invoke  basic  similarity  arguments. 
If  we  assume  that  the  potential  energy,  which  derives  from  the  wake 
kinetic  energy,  scales  by  the  external  parameters  U  and  H0,  then  we 
have  that 


j/S  "SV3 


so  that  the  scaling  requires 


F~' 

'  w 


4  4 


=  cons 


where  F„ 


is  the  inverse  internal  Froude  number  based  on  the  wake 


(33) 


characteristic  scale  H0.  Thus,  we  arrive  at  the  relationship 


ew  oc  (%)*/?  f^/n) 


for  the  normalized  rate  of  internal  wave  energy  production  by  the  wake 
collapse  mechanism. 

In  order  to  evaluate  the  .onstant  of  proportionality,  recourse 
must  be  maos  to  empirical  data.  Fig.  6  is  a  plot  of  fg  as  a  function 
of  Fw~ ^ ,  baoed  on  available  data  (Schooley  and  Stewart,  1963,  Stock¬ 
hausen  et  _al,  1966;  unpublished  data  provided  by  K.  G.  Williams  of  NRL) 
from  experimental  studies  with  self-propelled  bodies*,  the  self-propul¬ 
sion  mode  being  the  most  physically  relevant.  Although  there  is  no 
theoretical  justification  for  such  a  choice,  the  data  is  consistent 
with  a  relationship  of  the  form 


£  =  /.  S  K 


Examination  of  Fig.  7,  which  illustrates  normalized  vertical  density 
profiles  near  the  vertical  centerline  of  a  wake  at  maximum  vertical 
expansion  based  on  the  data  in  the  report  by  Stockhausen  jit  jil  (1966) , 
indicates  that  an  equivalent  uniformly  mixed  wake  would  have  Nj2  *1/2, 
for  Fw“l  96  0.04.  Thus,  based  on  the  relationships  (32)  and  (34)  and 
this  single  piece  of  data,  we  obtain  the  approximate  relationship 


=  o./  Fv 


It  is  important  to  emphasize  that  there  is  no  theoretical  basis  for  the 
expressions  (34)  and  (35),  although  they  are  consistent  with  the  avail¬ 
able  data  and  the  relationship  (32)  derived  from  similarity  arguments. 

In  deriving  (32) ,  what  in  fact  we  have  done  is  to  replace  the  generating 
body  with  a  virtual  wake  source,  so  that  when  results  which  are  incon¬ 
sistent  with  the  reality  of  the  body  are  encountered,  the  validity  of 
the  model  comes  into  question.  The  basic  constraints  which  must  be 
satisfied  if  the  model  is  to  be  at  least  internally  consistent  are  that 
fjj  <*  1  and  fg  ^  1,  which  introduces  the  bounds 


0.0/ 


is- 


Within  these  bounds,  we  then  have  an  approximate  expression 

£*  =  °-ooz 

}|c 

The  experiments  by  Williams  involved  simulated  self-propulsion,  in 
which  a  model  with  attached  propeller  was  towed  thru  a  tank,  the  pro¬ 
peller  being  driven  at  a  rate  such  that  the  pitch  of  the  propeller 
times  the  rotation  rate  was  equal  to  the  forward  speed  of  the  model. 
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for  the  rate  of  internal  wave  energy  production  by  the  wake  collapse 
process. 

The  normalized  rate  of  energy  supply  to  the  internal  wave  field 
as  determined  above  is  illustrated  in  Fig.  8.  Included  in  Fig.  8  is 
an  indication  of  the  relative  location  of  the  data  points  which  went 
into  determining  the  expressions  (34)  and  (35).  The  important  feature 
of  the  assumed  functional  form  for  the  relationship  (33)  is  the  inverse 
proportionality  of  Ey,  and  To  /  7^  ,  which  is  at  least  consistent  with 
the  idea  that  a  faster  body  will  create  a  more  energetic  wake,  and  the 
expression  (37)  should  provide  a  reasonable  estimate  of  the  rate  of 
internal  wave  generation,  at  least  over  the  parameter  ranges 
covered  by  the  model  studies. 


The  displacement  thickness  effect 


The  presence  of  the  boundary  layer  of  retarded  fluid  adjacent  to 
the  body  causes  an  outward  displacement  of  streamlines  in  the  surrounding 
fluid  which  is  perpetuated  by  the  turbulent  wake  and  results  in  an 
effective  increase  in  the  size  of  the  body  over  that  presented  to  the 
system  in  the  absence  of  viscous  effects.  The  dynamics  of  such  near 
field  flow  processes  will  be  largely  decoupled  from  the  ambient  strati¬ 
fication  provided  the  characteristic  time  scale  9*  of  the  process  is 
small  compared  to  that  associated  with  the  environmental  density 
gradient  7%  =  IT1.  If  we  assume  that  7$  «  7£,  and  if  the  body  is 
essentially  axisymmetric,  then  the  mean  flow  in  the  near  field  may  be 
described  by  the  Stokes  stream  function  ^  ,  where 

U  =  -  £  (38a) 

r  3  r 


V  - 


±  dJt 

K  3X 


(38b) 


where  u  and  v  are  the  axial  (x  direction)  and  radial  (r  direction,  where 
r2  =  y2  +  z2)  velocity  components.  Cylindrically  shaped  surfaces  of 
constant  are  the  Stokes  stream  surfaces,  which  are  the  axisymmetric 
analog  of  streamlines  in  two-dimensional  flows. 


A  number  of  parameters  have  been  developed  to  describe  the  proper¬ 
ties  of  boundary  layers  and  wakes.  For  the  purposes  of  this  report,  two 
integral  parameters  are  of  interest,  the  displacement  thickness  <St  ,  and 
the  momentum  loss  thickness  ,  where  we  define* 


*  These  definitions  differ  somewhat  from  the  more  traditional  defini¬ 
tions  which  are  carry-overs  from  two  dimensional  theory.  The  present 
definitions  are  used  to  facilitate  interpretation  in  the  sequel.  A 
thorough  discussion  of  the  concept  of  displacement  thickness  is  given 
byLighthill  (1958). 
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(39a) 


oo 

1  OS,"  =  fr(U-u)  dr 
*  'i 

00 

1  U  *  £  *  -jh  u  (O'-  u) 


(39b) 


where  rj,  is  the  local  body  radius  and  u  is  the  axial  flow  speed  in  the 
boundary  layer.  Since  we  are  primarily  interested  in  the  effects  asso¬ 
ciated  with  the  turbulent  boundary  layer  and  its  extension  as  the  turbu¬ 
lent  wake,  it  is  convenient  and  conceptually  advantageous  to  replace 
the  body  with  a  momentum  source  along  r  =  0  which  will  produce  an  equiva¬ 
lent  boundary.  In  this  case,  r^,  =  0  and,  to  the  boundary  layer  approxi¬ 
mation*,  the  displacement  thickness^  may  be  interpreted  as  an  effective 
outward  displacement  of  the  Stokes  stream  surface  =  0,  insofar  as  the 
flow  outside  the  boundary  layer  is  concerned.  To  illustrate  this  inter¬ 
pretation,  consider  a  flow  such  that  u  ••=  U  for  t  >& .  For  /-  >  S  ,  the 
stream  function  is  giwen  by  t- 

p  *  -  jr’u  dr'  -  Jr'  O  dr' 

-  kus,"- k  Uh* 

If,  however,  the  stream  surface  =  0  were  situated  at/-  =  <£ ,  instead 
of  r  =  O,  and  the  flow  ^vere  uniform  (£/«£/),  we  would  have 


O  dr‘ 


-  £  US,"  -  £  Or* 

Thus,  the  displacement  thickness  is  interpreted  as  the  radius  of  an 
equivalent  body  in  potential  flow.  This  is  the  "displacement  body". 

The  momentum  loss  thickness  is  a  measure  of  the  frictional  drag  on 
the  upstream  portion  of  the  body. 

A  number  of  techniques  for  computing  turbulent  boundary  layers  in 
general  are  available  (see,  for  example,  Kline  et  al,  1969)  some  of 
which  possess  fairly  straightforward  generalizations  to  boundary  layers 
on  axisymmetric  bodies  (see,  for  example,  Tsakonas  and  Jacobs,  1960, 
Yamaoka,  1963,  Nelson,  1964,  1966,  Tetervin,  1969),  If  the  body  is 
self-propelled,  the  computation  is  rendered  considerably  more  compli¬ 
cated  by  the  pressure  gradients  near  the  tail  of  the  body  induced  by 


See  sections  Ila  and  XlXa  of  Schlichting's  (1968)  bock  for  a  discussion 
of  the  boundary  layer  approximation. 
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the  propeller.  In  this  case,  iterative  procedures  (see,  for  example, 
Amtsberg,  1960,  Hucho,  1969;  see  also  Beveridge,  1969)  may  be  used  to 
compute  boundary  layers  which  compare  favorably  with  experiment  over 
all  but  the  very  tail  end  of  the  body,  in  the  immediate  vicinity  of 
the  propeller. 

Behind  the  body,  the  turbulent  boundary  layer  is  extended  as  a 
turbulent  wake.  In  the  absence  of  self-propulsion,  the  wake  is  charac¬ 
terized  to  the  boundary  layer  approximation  by  a  constant  momentum  loss 
thickness.  The  condition  of  self-propulsion,  on  the  other  hand,  dras¬ 
tically  alters  the  integral  properties  of  the  boundary  layer.  Again 
to  the  boundary  layer  approximation,  the  condition  of  self  propulsion 
requires  that  c^=  0,  that  is,  that  there  is  no  net  drag  on  the  body- 
propeller  system.  A  complete  description  of  the  interaction  between 
the  boundary  layer  and  the  propeller  is  not  possible  with  the  present 
state  of  the  art  of  hydrodynamic  theory.  However,  since  we  are  pri¬ 
marily  interested  in  the  displacement  thickness  of  the  wake,  which  is 
an  integral  property  and  hence  should  be  somewhat  insensitive  to  the 
details  of  the  flow  structure,  we  may  estimate  the  change  in  <5,  on 
passing  thru  the  region  of  strong  propeller  influence  in  a  fairly 
simple  manner,  using  the  condition  of  self  propulsion. 


A  slight  distance  upstream  of  the  propeller,  the  displacement 
and  momentum  thicknesses  are  S/  and  <fA  ,  while  downstream,  the 
displacement  thickness  is  ,  with  =  0.  Representing  the 

downstream  velocity  as  U  +  Up  ,  where  u  is  the  upstream  velocity 
profile  and  Up  is  the  modification  thereof  induced  by  the  propeller. 


we  have 


^  US*  =  / (U’  u)relr 

O 

£U*SAZ  Ju(U-u)rcJr 

o  90 

£  U  cf,p  -  / (U  -U- up)  h 
o  oo 


v 


If  we  represent  up  as  a  uniform  velocity  over  the  upstream  boundary', 
layer  thickness  £  ,  where  U  -  u  =  0  for  r  >  (5  ,  then  we  have  that  \ 

josS  -  i  os;  \ 

*{(/%** 
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The  condition  of  self  propulsion  (  S„  =  0)  then  requires  that 


(40) 

In  practice,  S  may  be  of  order  30%  to  50%  of  the  maximum  body 
radius  H,  6/  may  be  15%  to  25%  of  H,  and  the  shape  parameter 
Hr\  =  Si  /  6*  will  generally  be  greater  than  1.3  and  may  be  greater 
than  2  (Burstein,  1965,  Reed  et  _al,  1966,  Hucho,  1969),  so  that  up 
may  be  from  7%  to  10%  of  U,  and  S/a  may  be  from  10%  to  20%  of  the 
maximum  body  radius.* 

As  noted  previously,  the  dynamics  of  turbulent  wake  development 
in  a  stratified  medium  is  not  fully  understood.  However,  for  the 
momentumless  wake  in  a  homogeneous  fluid,  an  extensive  experimental 
investigation  exists,  due  to  Naudascher  (1965).  At  sufficiently 
large  distances  downstream  of  the  generating  element,  Naudascher 
found  an  essentially  self-preserving  wake  form.  Taking  advantage  of 
the  self-preserving  quality  of  the  wake,  we  may  use  Naudascher' s  data 
to  examine  the  change  in  displacement  thickness  with  distance  down¬ 
stream.  Fig.  9  illustrates  the  decay  of  the  normalized  displacement 
thickness  A  S, /H0  as  a  function  of  nondimensional  distance  downstream 
x/Hoj  where  A  is  a  constant  which  depends  on  the  exact  shape  of  the 
self-preserving  form  of  the  wake.  Extrapolation  of  Naudascher* s  data 
indicates  that  the  displacement  thickness  is  reduced  essentially  to 
zero  at  a  downstream  distance  x/Ho  =  150.  This  indicates  that  for  a 
self-propelled  body  with  length/diameter  aspect  ratio  R  =  7  1/2  and 
H0/H  =  1/2  in  a  homogeneous  fluid,  the  displacement  body  associated 
with  the  turbulent  wake  would  be  approximately  5  tines  as  long  as  the 
generating  body,  although  as  Indicated  previously,  the  maximum  radius 
of  the  displacement  body  would  be  only  10%  or  20%  of  the  maximum 
radius  of  the  generating  body. 

In  order  to  assess  the  relative  importance  of  the  displacement 
thickness  effect,  we  must  consider  its  interaction  with  the  hull  and 
wake  collapse  effects.  In  obtaining  estimates  of  the  size  of  the 
effective  displacement  body  associated  with  the  turbulent  boundary 
layer  and  wake,  we  have  been  forced  to  rely  on  data  and  arguments 
which  are  strictly  applicable  only  to  homogeneous  fluids.  If  the 
time  scale  7#  =  h*A J,  where  L*  is  the  half-length  of  the  effective 
displacement  body,  becomes  comparable  to  7%  ,  we  may  expect  some  dis¬ 
tortion  of  the  displacement  body  due  to  buoyancy  effects.  This  distor¬ 
tion  should  be  reflected  in  a  distortion  of  the  turbulent  wake  structure 
which  is  similar  to  but  distinct  from  that  associated  with  the  wake 

For  model  studies  in  which  the  propeller  radius  is  nearly  equal  to 
the  body  radius,  <£^raay  be  a  considerably  larger  fraction  of  the 
body  radius. 
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collapse  mechanism,  although  in  general  we  would  expect  the  effect  of 
this  interaction  to  be  of  minor  importance.  Potentially  more  important 
is  the  interaction  with  the  hull  effect.  An  examination  of  Fig.  4 
indicates  that  the  presence  of  an  effective  extension  of  the  generating 
body  will  enhance  the  supply  of  energy  to  the  hull  effect  internal  wave 
field.  The  exact  magnitude  of  the  relative  increase  in  hull  wave  making 
propensity  would  be  difficult  to  assess  with  the  relatively  crr.de  models 
discussed  in  this  report.  For  the  parameter  range  considered  above 
(maximum  displacement  thickness  of  approximately  1/2  H0,  H0/H  ==  1/4) 

the  increase  should  be  slight,  although  for  model  studies  for  which 
H0/H  is  generally  of  order  unity  the  effect  may  be  sensible,  and  signif¬ 
icant  increases  may  be  found  for  non-self  propelled  models  for  which, 
due  to  the  absence  of  the  propeller  contraction,  the  maximum  displace¬ 
ment  thickness  is  greater,  and  for  which  the  characteristic  length  of 
the  displacement  body  is  substantially  increased  (see,  for  example, 
Tennekes  and  Lumley,  1972,  sec.  4.2  and  4.3). 

Conclusions 

We  have  systematically  investigated  the  basic  physical  processes 
which  are  associated  with  the  passage  of  a  body  thru  a  stratified 
fluid  and  which  result  in  disturbance  motions  which  are  steady  in  a 
reference  frame  moving  with  the  body.  Three  dominant  mechanisms  are 
involved  in  the  flux  of  energy  into  the  internal  wave  field  which  is 
stationary  with  respect  to  the  body.  The  first  two  mechanisms  (the 
hull  effect  and  the  wake  collapse  effect)  act  somewhat  independently, 
and  have  different  parametric  spheres  of  influence.  For  a  given  body 
shape  and  size  and  a  given  environmental  stratification,  the  body 
effect  dominates  he  wave  field  at  relatively  low  speeds,  while  the 
wake  collapse  effect  dominates  at  relatively  high  speeds.  For  example, 
if  we  consider  a  body  with  character isitc  half-length  L  =  50  m  and 
length-breadth  aspect  ratio  R  =  10,  travelling  in  a  fluid  with 
N  «  10T2sec_1,  then,  assuming  a  propulsive  system  such  that  the 
characteristic  wake  radius  is  1/2  the  characteristic  body  radius,  we 
may  expect  the  body  effect  to  dominate  the  radiated  wave  field  for 
speeds  less  than  approximately  1  m/sec,  while  the  wake  collapse  effect 
should  dominate  for  speeds  in  excess  of  2  ra/sec.  The  third  effect  is 
important  insofar  as  it  interacts  with  the  other  effects,  the  dominant 
interaction  occurring  with  the  hull  effect.  This  mechanism,  the  dis¬ 
placement  thickness  effect,  produces  an  effective  extension  of  the 
body,  and  a  consequent  increase  in  the  supply  of  hull  effect  type 
energy  to  the  internal  wave  field.  Although  the  net  increase  in 
radioed  energy  due  to  this  effect  is  probably  of  minor  importance  in 
situations  similar  to  that  described  in  the  above  example,  it  may  be 
cf  significant  importance  in  model  studies  which  do  not  involve  self 
propulsion,  resulting  in  a  misleading  amplification  of  hull-type 
internal  wave  energy  in  such  situations. 


Appendix 


Problems  associated  with  the  bluff  sides  of  an  ellipsoid  may  be 
satisfactorily,  albeit  not  rigorously,  resolved  by  considering  a  new 
body  obtained  by  multiplying  the  local  height  of  the  ellipsoid  by  the 
vertical  projection  of  the  outward  normal  of  the  ellipsoid  at  that 
point.  The  intuitive  justification  of  such  an  approximation  arises 
from  the  fact  that  linear  theory  accounts  directly  for  only  the  verti¬ 
cal  component  of  induced  velocity.  Thus,  we  replace  the  body  shape 


by  the  new  shape 


A1 


for  which  the  induced  vertical  velocity  vanishes  around  the  edge. 


tfe  may  proceed  as  before, 
modified  body  shape 


h  (*J)  z  8rr  L 


obtaining  the  Fourier  transform  of  this 


c Z[L(x*  +  A*/R*)*J 

L*(X*+la/R*)  A2 


so  that  the  rate  of  internal  wave  energy  production  is 

/J  (xA  +  X*)v* 


O  O 


f  . / 


A3 


Nondimensionalizing  and  normalizing  as  before,  we  then  have  that 

£  z  *sj  /  ??  f/r' [(  £* /'  -  4V'* 

***  1  JJ 


A4 
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This  expression  would  then  replace  the  expression  (21)  for  the  normal¬ 
ized  rate  of  internal  wave  energy  production  associated  with  the  hull 
effect. 
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Tu/Tb 


Fig.  2  -  Normalized  internal  wave  energy  production  rates  (E)  as  a 
function  of  the  ratio  TJ,  /T^  =  NL/U  for  various  three 
dimensional  bodies. 
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Anomalous  behavior  of  normalized  energy  production  rate 
(E)  for  small  values  of  the  ratio  Tu/Tb  *>  NL/U  for 
ellipsoids  with  aspect  ratio  R.  It  is  expected  that 
production  rates  illustrated  in  this  figure  are  generally 
larger  than  those  actually  realized  (see  text). 
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Fig.  4  -  Normalized  internal  wave  energy  production  rate  (E)  as  a 
function  of  the  ratio  T^/  =  NL/U  for  ellipsoids 

with  aspect  ratio  R.  Production  rates  illustrated  in 
this  figure  are  expected  to  be  reliable  provided  that 

R_1TU/Tb  <  1/2. 
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Idealized  model  of  cross-section  of  partially  mixed 
wake,  as  defined  by  isopychnals,  illustrating  the 
displacement  function  6  .  Stratification  inside 
mixed  region  has  Brunt-Vaisala  frequency  N*,  that  in 
environment  has  N. 
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Wake  expansion  factor  (fE  =  Ae/Aq,  where  AE  is  the 
cross-sectional  area  at  maximum  vertical  expansion  and 
AQ  is  the  initial  or  characteristic  area)  as  a  function 
of  Ftf  =  NH0/U,  where  H0  is  the  initial  radius.  The 
dashed  line  has  a  slope  of  -1/2. 

0 “  Schooley  &  Stewart  (1963) 

A  -  Stokhausen,  Clark  &  Kennedy  (1966) 
o  -  Williams  (unpublished) 


Fig.  8  -  Normalized  power  Ofy)  associated  with  wake  potential 
energy  density  as  a  function  of =  NL/U,  for 
various  body  aspect  ratios  (R  =  L/H) ,  where  H/Hq  is 
the  ratio  of  body  radius  to  characteristic  initial 
wake  radius.  Normalization  is  the  same  as  that  used 
for  "hull  effect"  internal  wave  energy  production 
rate.  Points  indicate  parametric  location,  of  data 
used  in  evaluating  the  constant  of  proportionality 
in  expression  (31). 
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Normalized  displacement  thickness  A  5  ^/Ho  as  a  function  of  nondimensional 
distance  downstream  for  a  momentumless  wake,  where  A  is  a  constant  determined 
by  the  wake  shape,  based  on  data  from  Naudascher  (1965). 


